#!/bin/bash
#
# The purpose of this script is to find pairs of IRAC1 and IRAC2 files
# in a certain folder, find the corresponding CAT file in another
# folder, and run APT on each selected file.
# 
# Assumptions:
#  1. Your APT preferences file has the proper analytical parameters.
#  2. It's okay for the script to update the contents of the "photometryFile"
#     parameter in the APT preferences file.  (The script names the 
#     output .tbl file from APT according to the name of the input fits file.)
#
# The script quits if APT exits with a nonzero status.
#
# Author:
#  Mark Abajian               (626) 395-1811
#  Infrared Processing and Analysis Center
#  California Institute of Technology
#  Pasadena, California
#
# Date:  2011 August 16
#

script_basename=`basename $0`

function usage() {
  echo "Usage: $script_basename"
}

if [ $# -ge 1 ] ; then
    usage
    exit 1
fi

#for fitsfile in ../Data/NITARP-TRS1/*_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS1/3C175.1_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS1/3C[23]*_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS1/[^3]*_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS2/*_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS3/*_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS3/SDSS1{5[^0],6}*_irac[12].fits ; do
for fitsfile in ../Data/NITARP-TRS3/SDSS{150759,165943}*_irac[12].fits ; do
#for fitsfile in ../Data/NITARP-TRS4/*_irac[12].fits ; do
    photfile=`basename $fitsfile .fits`.tbl
    if [ `basename $fitsfile _irac1.fits` == `basename $fitsfile` ] ; then
	# the filename is not _irac1.fits, it must be _irac2.fits
	catfile=TRS168_irac12_cat/`basename $fitsfile _irac2.fits`_irac12.txt
    else
	# the filename is not _irac1.fits
	catfile=TRS168_irac12_cat/`basename $fitsfile _irac1.fits`_irac12.txt
    fi

    # Edit the preferences file to set the correct photometry file name
    preffile=$HOME/.AperturePhotometryTool/APT.pref
    sed -i -e "/photometryFile/s:/[^/]*$:/$photfile:" $preffile 

    echo "#########################################################################"
    echo "#########################################################################"
    echo "####    $fitsfile   ####"
    echo "#########################################################################"
    echo "#########################################################################"
    echo APT.csh $fitsfile $catfile
    APT.csh $fitsfile $catfile

    if [ $? -ne 0 ] ; then
	echo "Error:  Some problem running APT.csh $fitsfile $catfile"
	exit 1
    fi
done
